Kapitza homogenization of deep gratings for designing dielectric metamaterials 
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We theoretically investigate the homogenization of the dielectric response to transverse electric 
waves of a transverse grating characterized by the Kapitza condition, i.e. the permittivity is rapidly 
modulated with a modulation depth scaling as the large wavelength to modulation period ratio. 
We show that the resulting effective dielectric permittivity, in addition to the standard average of 
the underlying dielectric profile, has a further contribution arising from the fast and deep dielectric 
modulation. Such contribution turns out to be comparable with the other one and hence can provide 
an additional way for designing dielectric metamaterials. As an example, we discuss an effective 
metal-to-dielectric transition produced by the Kapitza contribution obtained by changing the grating 
depth. 
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Materials patterned with a periodic dielectric permit- 
tivity are a major platform for achieving electromagnetic 
radiation steering and manipulation If the modu- 
lation period is much smaller than the wavelength, the 
periodic medium displays an homogeneous dielectric re- 
sponse in turn described by effective material parameters 
0. In the metamaterial context, subwavelength layered 
metal-dielectric structures have attracted a great atten- 
tion due to the number of exotic electromagnetic phe- 
nomena they host resulting from the possibility of tai- 
loring both the magnitude and the sign of the effective 
permittivity. Leading examples of such unusual phenom- 
ena are subwavelength imaging 0, 0] , negative refraction 
[Ej], hyperbolic response (6|, optical non- locality 0, 0, 
transverse power flow reversing [9] and transmissivity di- 
rectional hysteresis [To| in extreme nonlinear metama- 
terials. Generally, the response of a dielectric stratified 
composite is theoretically described by the standard ef- 
fective medium theory [2] according to which the ordi- 
nary and extraordinary effective permittivities are the 
weighted arithmetic and the harmonic means of the un- 
derlying layer permittivities, respectively. On the other 
hand, in Ref.[7], Elser et al proposed an improvement 
of such standard effective medium approach by showing 
that the optical response of a realistic multi-layered nano- 
composite can be strongly affected by non-local effects. 
Recently, Rizza et al have suggested that a novel 
propagation regime, not encompassed by the standard ef- 
fective medium predictions, occurs if a Transverse Mag- 
netic (TM) wave propagates along the modulation direc- 
tion of a very deep and subwavelength grating. Here the 
large grating depth plays a fundamental role and, more 
precisely, if such depth scales as the large ratio between 
wavelength and the grating period (Kapitza condition), 
the authors show that a different mechanism support- 
ing diffractionless propagation occurs which is not signif- 
icantly hampered by medium losses. 

In this Letter, we theoretically investigate the propaga- 
tion of Transverse Electric (TE) waves in a realistic trans- 
verse dielectric grating fulfilling the Kapitza condition. 
We show that the rapid grating with a large depth entails 
a medium homogenization which can not be described 



by the standard effective medium theory. Accordingly, 
in addition to the average of the underlying permittivity 
profile, the ensuing dielectric effective permittivity has a 
further contribution which is generally not a small per- 
turbation. In addition, we study the transition from the 
standard effective medium regime to the Kapitza one by 
keeping fixed the grating period and gradually increasing 
the modulation depth and we show that, in the consid- 
ered example, the dielectric Kapitza correction turns the 
metallic medium behavior into a dielectric one. The con- 
sidered electromagnetic situation has a stringent analogy 
with the classical mechanical motion of a particle sub- 
jected both to a conservative field (described by a poten- 
tial energy U with a characteristic variation time T) and 
to a rapidly varying force (oscillating with a frequency 

n) mm. in this 

case, the condition QT >> 1 entails 
a separation between the slow and fast dynamical scales 
and, in complete analogy with the dielectric effective per- 
mittivity of the Kapitza medium, the effective potential 
energy driving the slow motion is the superposition of 
the average of U and a contribution resulting from the 
underlying fast dynamics. Such approach has been ex- 
tended in several physical context and in particular in 
optics for achieving nonlinearity management [3, [H| . 

Let us consider a monochromatic TE electromag- 
netic field E(x, z) = Re {E(x, z) exp [—it (2ttc) /A]} e y (of 
wavelength A) propagating in a periodic medium and sat- 
isfying the Helmholtz equation 



V 2 E + kle(x)E = 0, 



(1) 



where the period of the dielectric permittivity e(x) is A 
and fco = 27r/A (V = e x d x + e z d z ). Specifically, we 
consider a medium whose relative dielectric permittivity 
admits the Fourier series expansion 
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(2) 



and which satisfies the Kapitza condition, i.e. the spatial 
grating modulation is very fast (rj = A/ A <C 1) and the 
grating amplitude is very large, it scaling as the inverse of 



FIG. 1: (Color on-line) Sketch of a TE plane wave impinging 
onto a Kapitza slab of thickness L with an incidence angle 0. 
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77. Following the Kapitza approach [12] and exploiting the 
multiscale asymptotic expansion, we look for a solution 
of Eq. CQ) as 

E(x, X, z) = E(x, z) +rj ^ E n (x, z) exp (ink^X) , (3) 

where X = is a fast coordinate, whereas i5 and the 
last term in the RHS are the slowly and rapidly vary- 
ing parts of the electric field, respectively. Substituting 
the ansatz of Eq.((3j) into Eq.(pQ) and selecting the terms 
with the same spatial frequencies (hence dropping the 
fast scale), we obtain the coupled equations 

V 2 E + k 2 e m E + k%J2 ( b ~n + Va-n) E n = 0, 

(b n E - n 2 £ n ) + 0(77) = (4) 

where n / and 0(rj)/r) is a finite complex number in 
the limit 77 — ^ 0. In the asymptotical limit 77 — ^ 0, the 
second of Eqs.(j4]) turns into a simple algebraic equation, 
i.e. E n = bn/n 2 E so that the first of Eqs.((4]) yields, for 
rj -> 0, V 2 E + k 2 e ef ?E = 0, where 

Note that, in analogy with the mechanical situation, the 
obtained effective dielectric permittivity of Eq.(j5]), in ad- 
dition to the weighted arithmetic mean of the permit- 
tivity profile e m (as predicted by the standard effective 
medium theory) has a further contribution produced by 
the small (but fast) component of the electric field r\ E n , 
and this term can significantly affect the dynamics of the 
slowly average and leading order field E. 

In order to check the goodness of the multiscale asymp- 
totic expansion and to discuss the impact of the ac- 
tual small value of 77, we consider TE waves impinging 
onto a Kapitza slab of thickness L as reported in Fig.l. 
Specifically, as a theoretical benchmark, we assume the 
slab grating to be described by the permittivity 

e = e m + Se cos [(2tt/A)x] . (6) 



FIG. 2: Comparison between numerical slab transmissivity 
T and reflectivity R and corresponding Kapitza quantities 
T e ff and i? e ff- (a) normally imping waves for different grating 
depths; (b) inclined waves for Re(Se) = 10. 



We have numerically evaluated the transmissivity and 
the reflectivity of the slab by resorting to 2D full-wave 
simulations performed with the comsol RF module [l6j . 
In our simulations, we have considered two vacuum lay- 
ers placed at the facets of the slab for providing external 
excitation and hence for evaluating the scattering param- 
eters (transmissivity and reflectivity). We analyzed only 
one period along the x-axis imposing Floquet phase shift 
boundary condition on the sides of the unit cell. Fur- 
thermore, we have used matched boundary conditions at 
the entrance and exit facets (orthogonal to the z-axis) 
of the integration domain [16]. In Fig. 2, we report the 
numerical transmissivity T = |E t | 2 /|Ei| 2 and reflectiv- 
ity R = |E r | 2 /|Ei| 2 labelled with square and circle solid 
lines, respectively (see Fig.l for the definitions of the 
incident Ei, reflected E r , and transmitted E t field am- 
plitudes, respectively) with A = 100 /im, A = A/10, 
e m -0.25 + z0.0125, Im(Se) = 0.01 and L = 300 
/am. We plot T and R as functions of the grating am- 
plitude Re(5e) for normally impinging waves (panel (a)) 
and as functions of the incidence angle for Re(Se) — 10 
(panel (b)). On the other hand, using the results of the 
multiscale approach, we have analytically evaluated the 
transmissivity T e ff and the reflectivity R e ^ by regarding 
the slab as homogeneous with the effective permittiv- 
ity e e ff = e m + [r]Re (Se)] 2 /2 evaluated by using Eq.© 
for the profile of the dielectric constant of Eq.(j6]). For 
comparison purposes, in Fig. 2 (a) and Fig. 2(b), we also 
plot T e ff (dashed lines) and R e ^ (dot-dashed lines). It is 
evident that the agreement between the numerical sim- 
ulations and the predictions of the Kapitza approach is 
quite remarkable even for the chosen not so small ra- 
tio A/ A = 1/10. From Fig. 2(a) we note that, for 
small values of Re(5e), both T and R slightly depend 
on the grating depth as predicted by the standard ef- 
fective medium theory. On the other hand, a threshold 
exists (Re (Se) ~ 7) after which T and R show a marked 
dependence on Re (Se) proving that the Kapitza regime 
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FIG. 3: Comparison between numerical transmissivity (a) and 
reflectivity (b) (solid lines) and those predicted by Kapitza ap- 
proach (stars) as functions of and 1/77. The semitransparent 
surface interpolating the results of the Kapitza approach is re- 
ported for clarity purposes. 

has been reached. In addition, Fig. 2 (a) also shows that 
the Kapitza contribution to e e ff can generally play an im- 
portant role since, in the considered situation, by vary- 
ing the grating amplitude Re(5e), we observe a metal- 
to-dielectric transition: the real part of effective dielec- 
tric permittivity is negative and positive in the region 
< Re (Se) < 7 (shadow area) and in the region where 
Re(5e) > 7, respectively. Figure 2(b) shows that the 
Kapitza approach holds for all the incident angles. Note 
that, by comparing Eqs.(|6j) and (j2j) for the situation re- 
ported in Fig. 2(b), we have a±i = z0.005, b±\ = 0.5 



(a n = b n = for \n\ ^ 1) so that the Kapitza require- 
ment of the grating depth scaling as the grating period 
to wavelength ratio is numerically fulfilled (b±i being 
close to one). In Fig. 3, we plot the numerical trans- 
missivity (panel (a)) and reflectivity (panel (b)) along 
with those predicted by Kapitza approach letting both 
the grating depth and period to vary with the constraint 
Re(Se) = T]~ l = A/ A. Here the Kapitza requirement 
of the grating depth scaling as the grating period to 
wavelength ratio is evidently numerically fulfilled. Both 
T and R have been evaluated for different values of 77 
and 9 in the case A = 100 /im, e m = —0.25 + z0.0125, 
Im(Se) = 0.01, L = 300 /im. From Fig. 3, we note that 
the discrepancy between the numerical and Kapitza pre- 
dictions is significant for 77" 1 < 5 and it is negligible 
for increasing value of rf~ x , in agreement with the above 
discussed multiscale asymptotic analysis. 

In conclusion, we have studied the homogenization of 
the electromagnetic response of a transverse dielectric 
grating satisfying the Kapitza condition, i.e. the ratio 
between the wavelength and the grating period is much 
larger than one and the grating depth scales as this ra- 
tio. We have showed that the regime is described by a 
novel effective medium theory predicting for the effective 
permittivity a Kapitza contribution to the average per- 
mittivity which is generally not a small perturbation. We 
expect our approach to offer new avenues for the tailoring 
of metamaterial properties and specifically for the design 
of epsilon- near-zero metamaterials |l7j . 
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